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Abstract: We have explored the question of whether the renormahzation group running 
of the neutrino mixing parameters in the Minimal Supersymmetric Standard Model is 
detectable with ultra- high energy neutrinos from active galactic nuclei (AGN). We use as 
observables the ratios of neutrino fluxes produced at the AGN, focusing on four different 
neutrino production models: [^l^+p^ : ^>°^+p^ : = (1 : 2 : 0), (0 : 1 : 0), (1 : : 0), 

and (1:1:0). The prospects for observing deviations experimentally are taken into 
consideration, and we find out that it is necessary to impose a cut-off on the transferred 
momentum of > 10^ GeV^. However, this condition, together with the expected low 
value of the diffuse AGN neutrino flux, yields a negligible event rate at a km-scale Gerenkov 
detector such as IceGube. 
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1 Introduction 



Experiments performed over the last ten years have confirmed the mass-induced neutrino 
oscillation phenomenon, caused by the non-coincidence of mass and flavour eigenstates. 
Rather, the flavour states, z^e, i^^, and Vr, are linear superpositions of the mass eigenstates, 
which are customarily denoted by vi, V2, and z/3, i.e., 



(1.1) 



with a = e, fj,,T and i = 1,2,3. The coefficients Uai are components of the lepton mixing 
matrix, U,^, also known as the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, that 
connects the bases of mass and flavour neutrinos and which can be parametrised in terms 
of three mixing angles, 9i2, 023 and ^13, and one CP-violating phase, 6, as 
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where Cij = cos (Oij) and Sij = sin {6ij). A recent global analysis [1] used data from 
solar (SNO, Chlorine, Gallex/GNO, Borexino, SAGE), reactor (KamLAND, CHOOZ), 
accelerator (K2K, MINOS), and atmospheric (Super-Kamiokande) neutrino experiments 
to find the following 3a bounds on the mixing angles and mass differences: 



Amsi - 7.59_o.69 X 10 eV , ^"^3i - | ^^2.46 ± 0.37) x 10"^ eV 

+3.2° fl__/n 8+10.7° .__^„+6.9 



(1.3) 



0i2 = 34Atti , 023 = 42.8 ^^^3^ , 0i3 = 5.6t^:r, 6 e [0,27t] . (1.4) 

Note, however, that the experiments performed so far have probed neutrino energies from 
the few MeV (solar neutrinos) to the GeV range (atmospheric neutrinos). Thus, it is 
important to point out that the values of the mixing angles and mass differences found in 
this and similar analyses are valid, in principle, only in that energy range. 

Ever since the oscillation phenomenon was first observed, neutrino physics has been 
considered a window towards physics beyond the Standard Model (SM) [2]. Thus, the 
study of the influence of new physics on neutrino properties is of high importance. In the 
following, we shall address how loop corrections on neutrino interactions are affected by 
new physics, in particular that of the Minimal Supersymmetric Standard Model (MSSM). 

In neutrino experiments, loop corrections can play a role at two sources. The first 
source comes from the neutrino interaction vertex, where diagrams with intermediate par- 
ticles can modify the vertex by either introducing a correction due to form factors [3] or 
by changing the flavour structure of the mixing matrix entering the vertex itself (an ex- 
ample for the CKM matrix in the MSSM can be found in [4]). At one loop, most of the 
corrections in the SM are proportional to the PMNS matrix, so we do not expect to obtain 
any signiflcant corrections to the flavour structure. On the other hand, supersymmetric 
models like the MSSM introduce chargino, neutralino and slepton mixing matrices into the 
game, meaning that they have the potential to modify the flavour structure of the effective 
vertex. 

As the vertex corrections need to be renormalized, an arbitrary scale dependence is 
introduced. This causes the well-known running of parameters through renormalization 
group equations (RGEs). In order to avoid large logarithms in the loop functions, it is 
customary to set the scale /x such that these large logarithms vanish. For vertex corrections, 
this means that the scale is of the order of the transferred momentum, ~ Q = \/—q' [5]. 

The evaluation of the Green function related to the vertex requires all parameters 
involved in the function to be set at the same scale. This implies that one needs to take 
into account the RGE evolution of all parameters according to the transferred momentum 
of the interaction. Moreover, as the RGE evolution of the parameters is independent of 
the process one is analysing, this constitutes an independent second source of corrections 
to the vertex. The running of the mixing angles has been thoroughly studied, for example, 
by [6-10]. The difference in the evolution of the mixing parameters between the SM and 
MSSM was analysed in detail in [9, 10], where it was shown that if the bounds on neutrino 
masses are respected, the SM running is negligible. Thus, even though in this work we 
shall perform the SM running, we shall refer to it as the "no-running" scenario. The MSSM 
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scenario, on the other hand, could provide an additional enhancement that allowed large 
running effects to take place. 

Motivated by the results in [9], in this work we study whether there could be any 
observable deviations from the SM neutrino oscillation framework probed at a very high 
energy scale, due only to the RGEs in the MSSM. To this end, in Section 2 we briefly 
review the origin and structure of the RGE corrections. In Section 3 we establish how 
these corrections should be introduced within a neutrino experiment. In Section 4 we 
describe the changes one could expect if the RGE corrections are implemented from ultra 
high energy neutrinos produced by active galactic nuclei and examine the experimental 
requirements in order to observe the deviations. We conclude in Section 5. 

2 Running of Neutrino Mixing Parameters 

In a SM with massive neutrinos, the RGE evolution of neutrino masses and mixing angles 
is a consequence of the renormalization of the effective neutrino dimension-five operator 



where h.y is the scale where the new physics generating this operator is decoupled. When 
the electroweak symmetry is broken, this leads to the following mass matrix: 



where v is the vacuum expectation value of the Higgs field H. 

In SUSY models, we can build a SUSY operator analogous to Eq. (2.1) by replacing 
all fermion fields by superfields, and the H field by the SUSY Hu superfield. This modifies 
the mass matrix of Eq. (2.2) by a factor sin^ /3. Note that this relies on the assumption that 
either K^, is larger than the SUSY breaking scale, such that the corresponding dimension- 
five operator must be holomorphic, or that the non-holomorphic SUSY-breaking corrections 
are small. 

Although at high scales the neutrino mass matrix is highly model-dependent (i.e., see- 
saw mechanisms [11-14], radiative mass generation [15, 16], SUSY with R-Parity break- 
ing [17]), at lower scales the situation is different. After assuming a value for the lightest 
neutrino mass and the scale A^, one can uniquely define the dimension- five operator using 
only the neutrino masses and mixing parameters. This implies that, to study the RGE 
evolution of the neutrino parameters within models like the SM and MSSM, one can intro- 
duce the neutrino mass matrix operator at the electroweak scale and not worry about its 
origin until the moment one reaches the scale Ky. 

Furthermore, this operator is the only dimension-five operator allowed by the gauge 
symmetries. Given the dimensionality and flavour structure of this operator, its renormal- 
ization must always be proportional to the operator itself, leading to the following RGE 
(3u function [9]: 
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where x = ln(/////o)) A* is the scale where the operator is being evaluated, and fiQ is the 
scale where the initial conditions are specified. Here, and in the following equations, Yj 
shall denote the Yukawa matrix of the fermion f = u,d, e. The coefficients C and a were 
calculated in [6-8]: 

C=/-^ (2 4) 

\ 1 (MSSM) ^ • ^ 

^ ^ f -35i + 2 TviYhe) + 6 Tr(yjy„) + 6 TviY^Ya) + A (SM) 

\ -yi - Qgl + 6Tr(yJy,) (mssm) ^ ' ' 

A study of the variation of the neutrino mixing parameters in the SM and MSSM was 
carried out in [9, 10]. The authors found that, below the scale, the variation of the 
mixing angles could be approximated using the following analytical equations: 

til o I mi e*''^^ + 7712 6**^2 1 ^ 
012 = Sin 2^12 sin2 023^ L + 0(0^3) 

2 / _\_ Ynnl^ _ \m-i p'^^^ ' ''^ 



923 = -C:p;^ sin 2^23 cos^ 612 ' . 2 + ^12 \ 2 

?>1tx^ \ Amir, Ami 




= ^32^^ '^^^ '^^^ i - '^^ - A;;^ - '^^ - A777i,A7ni, 

+0(%) (2.6c) 

where y,- = (1^)33 and 777 j are the neutrino masses. The full set of analytical equations 
can be found in [9]. Moreover, although these equations are useful for understanding the 
variations obtained by the running, in all our results we shall solve the exact equations 
using the REAP package [10]. 

As these three equations are mainly dominated by y^, in the MSSM we can expect 
an enhancement of order {v/vdf' = (1 + tan^ /3) with respect to the SM. Thus, for values 
of tan ;S ~ 50, we can get very large deviations in the mixing parameters at high scales. 
Moreover, in all models, the derivatives of the mixing angles are proportional to the neutrino 
masses, so the variations can be further enhanced if the neutrino masses are large. A further 
source of enhancement comes from the Majorana phases, which can conspire in order to 
increase the effects of the running even further. 

As we want to investigate the largest possible variation of the neutrino flavour-transition 
probability with respect to its value evaluated at the best fit values, we shall maximize the 
777j/ and tan /? enhancements. In this work we shall use a normal mass hierarchy, and set the 
lightest neutrino mass at 7771 = 0.43 eV, which is the upper limit allowed by the WMAP- 
only 7- year data [18]. In addition, in the MSSM, we shall also use tan /3 = 50, and a SUSY 
decoupling scale AgusY = 1 TeV. 

3 Neutrino oscillation probability and RGE effects 

The analytical expressions shown in eq. (2.6) are of considerable use for model builders, as 
they can connect low-energy observables in the neutrino sector with high-scale theories like 
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Leptogenesis [19, 20], Flavour Symmetries [21-23], and Grand Unified Theories [24, 25]. 
However, it is not clear if the effects of the running could be observed in an experimental 
situation or not. 

As previously mentioned, when introducing quantum corrections to interaction ver- 
tices, the scale ^ is set equal to the transferred momentum Q = y—q^- On the other 
hand, deviations from the no-running scenario occur at scales larger than the SUSY de- 
coupling scale. Thus, we need large values of Q. These can be obtained, among other 
conditions, when we have attained large values of the neutrino energy, Ey. In particular, if 
we assume AsusY = 1 TeV, to get a Q of this order of magnitude we shall require energies 
much larger than those found in present and future long-baseline neutrino experiments. 
For instance, the largest energy considered for neutrino factories is of about 50 GeV. Thus, 
in order to reach larger energies, we must use astrophysical neutrinos. Examples of sources 
of these types of neutrinos are supernova remnants [26], gamma ray bursts [27], and active 
galactic nuclei ( AGN) [28] , with expected neutrino energies of order 10^^ 10^^ and 10^^ eV, 
respectively. For definiteness, we restrict our analysis to AGN, which provide the neutrinos 
with the highest possible energy in the Universe. 

In a neutrino experiment, there are two different moments where the scale is relevant: 
in the production and the detection of the neutrinos. In AGN, the neutrino production is 
realised though the decay of a particle, typically the pion, which sets a low scale for the 
transferred momentum (i.e., Q = m^). Meanwhile, when the ultra high energy neutrinos 
reaching the Earth interact through deep inelastic scattering with the nucleons of the 
detector, there is a possibility to achieve a large value of transferred momentum. It is in 
this last context where we can meet the necessary conditions to test the RGE effects in the 
neutrino system. Therefore, and given that we have two different Q scales, the Uy matrix 
used to calculate the incoming neutrino flux will differ from the Ul matrix that appears 
in the detection cross-section. We parametrise the Ul matrix, as in eq. (1.2), by three Cij 
mixing angles and a phase 5i, which are calculated by solving the RGE equations for the 
Oij mixing angles measured in the oscillation experiments. Due to the low scale involved 
in the neutrino production, the matrix is given by the current measurements of the 
neutrino oscillation parameters. These are shown in eq. (1.4), and from now on shall be 
those used in our analysis. 

As the Uy and Ul matrices can be factorised out of the production decay rate and 
detection cross section, we will be able to define a function Pq,/?, which will measure the 
probability of a neutrino produced with flavour a interacting with flavour /5: 



3 



(Q) = ^ !([/.)« [KiQ)) 
=1 



I Pi 



(3.1) 



Notice that if we set Ul ^ U^, this expression coincides with the standard expression for 
neutrino oscillation probability at large distances. 

It is useful to have an unraveled expression for this new probability. In the following, 
for conciseness, we shall write Qai = \{Ui,)ai\^- Making the dependence on the dj angles 



- 5 - 




Figure 1. Flavour-transition probabilities Pe/9, P^.^ versus difference between the Majorana phases, 
4>i — 4>2- Lines correspond to standard probabilities, with /3 = e (solid), fi (dashed), and r (dotted). 
Areas correspond to modified probabilities, where, at each point, (f>i has been randomly varied, 
while 02 has been adjusted to keep the desired phase difference. We have (3 = e (grey), fi (brown), 
and T (orange). The RGE evolution was performed at = 10^^ GcV^, fixing the dij and Amfj at 
their best-fit values, and (5 = 0. 



and the CP phase 5i explicit, we can write the probabilities as: 

Pae = QaS + C^ia [{gal " Oas) " S^^^ (Oal - ^02)] , (3.2) 
Pafi = 9a2 + {gal — 9a2) 

, f gal - ga2\ o\ 
+ ( 2 ] ^5i^2Ci2'52C23'5Ci3 ' (3-3) 

Par = ga2 + S^^j {^al — 5o2) 

+ ^23 [{9a3 - ga2) - S^^^ (l + S^w) {dal " S'a2) + 5^13 {9al " gas)] 

C5iS2Ci2'52C23'SCi3 ' (3-4) 

with c-^ = cosx and = sinx- The explicit dependence of the gai on the standard 
mixing angles can be found in Appendix A. Hereafter, we will consider only oscillation 
probabilities of the form P^p and {(3 = e,/Li, r), since the electron and muon neutrinos 
are the only ones that intervene in the different models of the neutrino production at AGN. 

Naturally, the probabilities do not have an explicit dependence on the Majorana phases. 
However, as the final value of the Qj angles depends on the value of each phase, we find 
that there does exist an implicit dependence on them. In particular, we have found that the 
higher the value of the Majorana phase difference 4'i~ 4'2, the higher the possible values of 
C12 and C13 that are respectively achievable. In ^23 there is a different behaviour: as (pi — </>2 
rises, the accessible region of C23 becomes focalised around a certain value, dependent on 
the chosen values for the mixing parameters at production. Note that, as long as 5 = 0, 
the regions are the same for (pi — (p2 and for (^2 — Nevertheless, using 5 7^ does not 
generate large deviations, and these statements remain valid. 

We find it illustrative to show how the probabilities depend on the Majorana phase 
difference, which is shown in figure 1. Here, the left and right panels show, respectively, the 
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Figure 2. Flavour-transition probabilities under MSSM running, as functions of the transferred 
momentum Q^. In the gray areas we vary all mixing parameters at production within their current 
3cr bounds (the CP and Majorana phases were varied between and 27r), while in the hatched areas 
only and the phases are varied. The dashed (red) bands indicate the limits corresponding to 
standard oscillations. 



probabilities Pe/3 and P^/j, for /3 = e (solid line, grey areas), fj, (dashed line, brown areas), 
and r (dotted line, orange areas). The lines correspond to the standard, no-running, 
probabilities, while the areas correspond to the probabilities influenced by the MSSM 
running. To obtain the areas, we have fixed all parameters at their best-fit values, and set 
5 = 0, = 10^^ GeV^. For each phase difference, we have varied between and 27r, 
and adjusted 02 to keep the difference at the desired value. 

Figure 1 shows a strong dependence of the transition probability on the Majorana 
phases, in clear distinction with the standard neutrino oscillation scenario. In fact, after 
the masses, the Majorana phases are the most important parameters in the neutrino mass 
operator affecting the final value of the probability. Note that the maximum deviation of 
the modified Pee and Per occurs at </>i — 02 = 0, while the maximum deviation of Pe^i occurs 
at </>i — 02 = ±71". The deviations are less pronounced for the probabilities, with P^e 
and P^r reaching maximum deviation at 0i — 02 = 0, and P^^ at 0i — 02 = ±71". 

It is also instructive to understand how the transition probabilities depend on the scale 
fi = Q, which is shown in figure 2. In this plot, the gray areas represent the accessible values 
of the probability after MSSM running when the standard mixing angles and squared- 
mass differences are varied within their current 3a experimental bounds, and the phases 



are varied between and 2tt. The hatched regions show the same, but considering only 
the variation of ^13 and the phases. From this we can understand the role played by the 
uncertainty in the other better-known parameters. The bands limited by dashed (red) lines 
correspond to the probability calculated using only standard oscillations, without MSSM 
running, when the mixing parameters are allowed to vary in the same manner. Notice that 
the deviations in the probabilities start at = 10^ GeV^ = Ag^gy and increase with Q^. 

Figure 2 allows us to make definite predictions regarding each probability. Although 
the high-scale behaviour of Pee is uncertain, we can expect Pefj, to generally increase, and 
Per to always decrease with Q^. In fact, we can obtain a suppression of Per down to 
zero, which shall provide very interesting results in the next sections. This behaviour is 
consistent with a large decrease of ^23) with deviations from the general trend depending 
on the value of the Major ana phases, as shown in figure 1. 

On the other hand, P^,- tends to increase with Q^. P^g and shall generally be 
expected to decrease, although it is possible for them to remain invariant. This again is 
consistent with a decrement of and deviations depending on the exact value of the 
Major ana phases. 

4 Astrophysical observables 

4.1 The UHE astrophysical neutrino flux 

AGN have long been presumed to be sites of ultra-high-energy (UHE) neutrino production. 
In the scenario of neutrino production by meson decay, it is assumed that within the AGN 
protons are accelerated through first-order Fermi shock acceleration [29] and that pious are 
produced in the processes 



The neutral pious decay into gamma rays through vr — ?■ 77, while the charged pions decay 
into electron- and muon-neutrinos through 

7r~^ ffj. + fj,'^ i^^, + e'^ + i^e , tt" V ^ + fj," + e~ + Ue + 1^^ . (4.2) 

If neutrinos are produced by these processes, then the ratios of the different flavours {ux + 
Vx) to the total fiux are: 



Note, however, that these flavour fluxes are approximate; a detailed analysis of the pion 
decay and the contribution from higher-energy processes will nevertheless result in values 
that are close to these standard ratios [30, 31]. 

By the time neutrinos reach Earth, neutrino oscillations will have redistributed the 
flux among the three flavours, so that the flavour fluxes at detection are 







(4.3) 




l3=e,fi,T 



(4.4) 
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For instance, standard mass-driven neutrino oscillations, evaluated at the best-fit values of 
the mixing parameters, eq. (1.4), distribute the total flux in the vr^ decay scenario in an 
approximately uniform manner among the three flavours, i.e., (^^p^+p^ : <I>i,^+p^ : ^ 
(0.36 : 0.33 : 0.31). 

Deviations from these values could signal the presence of new physics at work. In fact, 
there are several theoretical hypotheses, such as violation of Lorentz or CPT invariance [32, 
33] , that could in principle work on top of the standard mass-driven oscillation mechanism 
and induce these deviations. In particular, as we have shown before, the MSSM introduces 
RGE-induced changes in the standard neutrino oscillation probability, which implies that 
such deviations could be also achieved within this scenario. However, the expected flavour 
fluxes could be also modified by using a different neutrino production mechanism. Thus, 
it is essential to thoroughly describe the deviations we would expect from the MSSM, 
including simultaneously the different potential production processes we know of. Note 
that, while it is possible for the flavour ratios at production to have an energy dependence 
[34], we have not considered this possibility in our analysis. 

In a related production process [30, 34-36], the muons produced by pion decay may 
lose most of their energy before decaying, so that a puie-P^ flux is generated at the source, 
i.e., 

■■ : =(0:1:0). (4.5) 

Under standard oscillations, and with the mixing parameters set at their best-fit values, 
these initial fluxes yield, at Earth, {^u^+ue ■ ^u^+pp, '■ ^i^^+p^) ~ (0-26 : 0.36 : 0.38). 
Alternatively, a puie-Ue initial flux, i.e., 

■■ ■■ <-^Pr) =(1:0:0), (4.6) 

has been considered, e.g., in [30, 34]. In this scenario, high-energy nuclei emitted by the 
source have sufficient energy for photodisintegration to occur, but not enough to reach 
the threshold for pion photoproduction. The neutrons created in the process generate 
z7e through beta decay. Using standard oscillations and best-fit values for the mixing 
parameters, this yields, at Earth, ($p^ : : (0.55 : 0.26 : 0.19). 

Finally, semileptonic decays of charm quarks can generate flavour fluxes at production 
of [37, 38] 

: <I>°,+p, : '^l^p^) = (1/2 : 1/2 : 0) . (4.7) 

These fluxes can also be produced as the result of a pile-up effect [34] . In this scenario, stan- 
dard oscillations at best-flt values result in (<&,y^+p^ : <&i/^+p^ : ^u^+Pr) ^ (0.41 : 0.31 : 0.28). 

The initial neutrino and anti-neutrino flavour fluxes for the four production scenarios 
that we have considered are summarised in table 1. The fluxes at Earth for the four 
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Scenario 


Description 


($0 . 

\ l^e 


$0 . ^0 \ 

I^IJ. ■ i^T J 




I'll ■ i't y 


A 


Beta decay 


(0 


: : 0) 


(1 


: : 0) 


B 


Muon-damped 


(0: 


1/2 : 0) 


(0: 


1/2 : 0) 


C 


Charm decay 


(1/4 


: 1/4 : 0) 


(1/4 


: 1/4 : 0) 


D 


Pion decay 


(1/6 


: 1/3 : 0) 


(1/6 


: 1/3 : 0) 



Table 1. The four scenarios of neutrino flavour ratios at production time. The normalised fluxes 
<i>^^ and represent, respectively, tlie fraction of neutrinos and anti-neutrinos of flavour a pro- 
duced at the source. Note that, except for scenario A, the flux is equally divided between particles 
and anti-particles. 



scenarios are hence: 



Scenario A 



Scenario B : (E^) = -P^„ {E„) 



Scenario C : (E,) = - {P.^ (E,) + (E,)) 
Scenario D : <^,_. (E^) = ^ {Pea {E„) + 2P^„ (E^)) 



(4.8) 

(4.9) 
(4.10) 
(4.11) 



4.2 Flavour ratios 

Although we have four well-motivated models giving the flavour composition of the initial 
neutrino flux coming from AGN, the total neutrino flux is subject to many uncertainties 
(see [39] and references within). A good way to avoid having to deal with the uncertainties 
in the initial neutrino flux is using ratios of fluxes measured on Earth [40]. For instance, a 
typical ratio used is 



(4.12) 



which is related to the experimental ratio of Cerenkov events to the total number 

of neutrino events. If it is possible to distinguish between and v^- showers, another 
convenient ratio is 

R = ^J±±Ei. . (4.13) 



In what follows, we will show our expectations for R and T considering the four 
different scenarios, A, B, C and D, and assuming no running and MSSM running of the 
mixing parameters. 

In all these results, we variate ^13 in its 3a range, and all phases between and 27r. 
We will either hold all other oscillation parameters at their best-fit values, or variate them 
in their 3(7 range. Also, throughout the present section we shall fix = 10^^ GeV^ in 
order to maximise the MSSM RGE effects. 
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Figure 3. Expected ratio of fluxes assuming no running (orange) and MSSM RGE running (blue) 
of mixing parameters at ~ 10^^ GcV^, for a (1 : : 0) production model. We show only the 
variation of ^13 and the phases on the left, and the variation of all oscillation parameters on the 
right, as described in the text. 

4.2.1 Scenario A: Production through beta decay 

For simplicity, we start with the initial (1:0:0) ratio, which occurs if the neutrino 
production at AGN is governed by beta decay. We then have: 

r = PeM , ■ (4.14) 

From figure 2, we would generally expect an enhancement in Pg^, such that T is larger 
than the no-running typical values. Also, although whether Pee is enhanced or suppressed 
depends on the values of the oscillation parameters at production time, we have seen that 
Per is generally suppressed, with the possibility of being zero. Thus, we expect R to be 
able to reach very large values. 

In figure 3 we show the scatter plot of R vs. T . There we observe that the no-running 
expectations at Sa are T G (0.14,0.33) and R G (1,4), while the MSSM running expands 
this into T G (0.15, 0.47) and R G (2, 70). We see that it is possible for R to reach up to 18 
times the highest possible no-running value, a clear indication of new physics. T can also 
become very large, being particularly useful if the errors in other oscillation parameters are 
diminished. Note that the ranges of R and T can be extracted in a straightforward manner 
due to the directly or inversely proportional relationship between these observables and the 
oscillation probabilities. For instance, the upper (lower) limit of R has been obtained from 
the ratio of the maximum and minimum values of the oscillation probabilities in figure 2, 
PeelPe^i ~ 0.72/0.01 (0.4/0.2), at Q2 = iQii GeV2. 

It is interesting to note that since the probabilities shown in figure 2 are moderately 
fiat around = 10^^ GeV^, we can expect these results to hold even for somewhat lower 
values of Q^. This particular observation applies to the other three production scenarios 
as well. 

It can also be seen that the correlation of T and R is invaluable. Apart from being 
indicators of new physics, both ratios define a zone in which the new physics effect can be 
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identified as coming from SUSY. If either R or T are observed out of both orange and blue 
regions, the new physics origin of the deviations could not be attributed to MSSM RGEs. 

4.2.2 Scenario B: Muon-damped production 

The muon-damped production process generates the (0:1:0) fluxes. The T and R ratios 
can then be written as: 

T = P^^ , R=-^. (4.15) 

We can again use figure 2 to predict how these ratios shall behave. In contrast to 
the previous production process, we could expect some suppression in P^^, giving a value 
of T lower than the typical ones. Furthermore, the possible suppression of P^e and the 
enhancement of P^,- would also decrease R, giving a situation opposite to the (1:0:0) 
ratio. 

Figure 4 illustrates the situation. The no-running expectations at 3a are T G (0.33, 0.5) 
and R € (0.4,1.0), while the MSSM running takes this into T G (0.25,0.5) and R € 
(0.1, 1.0). As in the previous case, these ranges can be read directly from the corresponding 
range of P^^ , for T, and the proper combination of the maximum and minimum values of 
P^^ and Pfj,r, for R. 

In this case, the suppression of T and R is not as strong as their enhancement in 
the previous scenario. However, when correlated, it is clear that they can be disentangled 
from the no-running regions. Furthermore, it is also clear that the zone delimited by their 
correlation can be separated from the zone delimited in Scenario A. It should then be 
possible to deduce the original neutrino production process at the AGN. 

4.2.3 Scenario C: Production through charm decay 

The decay of D mesons can lead to a (1/2 : 1/2 : 0) production ratio. The detected ratios 
are then described by: 

T=liPe, + P,,) , R=^^^- (4.16) 
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Figure 5. Same as figure 3, but for a (1/2 : 1/2 : 0) production model. 



The results for T and R can thus be understood as an average of the predictions for 
muon-damped and beta decay production. As we have seen, these two scenarios have 
opposite behaviours, so it is possible for the deviations to cancel each other. In fact, 
for T, the possible reduction of P^^ and increase of Pg^ cause the total prediction to 
smooth out, getting milder maximum variations. For R, although Pgr is decreased, the 
large P^r does not allow R to reach the high values found in the beta decay production 
model. Furthermore, the smaller P^g combines with the somewhat larger P^e, such that 
the reduction characteristic of the muon-damped production model is not observed. All in 
all, R does not differ much from the no-running prediction. 

We show the no-running and MSSM predictions in figure 5. The no-running expecta- 
tions at 3(7 are T £ (0.30, 0.34) and R G (1.0, 1.8), while the MSSM running expands this 
into T e (0.26, 0.39) and R G (1.0, 2.5). We can see it is unlikely for neither T nor R to be 
used to distinguish new physics from no-running behaviour. By reducing the uncertainty 
in the mixing parameters one could hope to improve the situation, as the no-running region 
becomes very tiny, but, still, the difference between the no-running and MSSM running 
regions will not be as large as for scenario A. 

4.2.4 Scenario D: Production through pion decay 

The (1/3:2/3:0) model associated to full pion decay is arguably the one most studied in 
the literature. The ratios detected at Earth are described using: 

T=UPe, + 2P,,) , R=^^^^. (4.17) 

This is similar to the charm decay model, but since more weight is given to Pfifi, P^e 
and P^T, we can expect our results to lean slightly towards those for the muon-damped 
production model. Nonetheless, the variations are still not as significant as in the muon- 
damped and beta decay production models. 

This is what we observe in figure 6. The 3a no-running expectations are T G (0.32, 0.38) 
and R G (0.8,1.4), while the MSSM running expands this into T G (.27, .43) and R G 
(0.7, 1.7). No large deviations from the no-running expectations can be observed in any 
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Figure 6. Same as figure 3, but for a (1/3 : 2/3 : 0) produetion model. 



parameter, even if there is an improvement in the measurement of oscillation parameters, 
unless a large resolution in both T and R is achievable experimentally. 

4.2.5 Other Production Scenarios 

In the following we shall parametrise the initial ^^^g+p^ : ^Ufj,+u^ '■ ^ut+Pt^ as (1 : n : 
0), as in [38]. We show T and R as functions of n in figure 7, where we vary n from 10~^ to 
10^ in order to properly reproduce the limiting cases (1:0:0) and (0:1:0). Again, the 
orange region shows the no-running prediction for T and R, while the blue region shows 
that for the MSSM RGEs. 

From figure 7 we see that the previously analysed scenarios are representative cases 
when varying n, since the (1 : 2 : 0), (0 : 1 : 0), (1 : : 0), and (1:1:0) initial ratios 
correspond to n = 2, n = oo, n = and n = 1, respectively. 

However, figure 7 allows us to understand additional features of our framework. For 
instance, we can see that the no-running scenario is bounded. In particular, taking into 
account the 3a variation of all mixing parameters, T is roughly bounded within (0.15, 0.5) 
and R within (0.4,4). Any measurement of these parameters outside these bounds is a 
clear indication of new physics. 

Moreover, it is straightforward to see that the MSSM and no-running predictions for T 
cannot be disentangled if n is not known. A large deviation from the no-running scenario 
at ?i ^ is degenerate to a situation with no deviation at ?i — ?> oo, and vice- versa. Thus, 
on its own, T cannot provide any useful information for our purpose. 

In contrast, R unambiguously shows very large deviations from the no-running scenario 
whenever the condition | log]^o("')l ^ 1 is satisfied. In addition, the measurement of R not 
only can reveal the presence of new physics in the running of neutrino mixing angles, but 
it can also indicate the initial composition of the neutrino flux coming from AGN. For 
instance, if we measure R = 10, we not only know that there is new physics at work, but 
we can also set a limit of n < 0.2. Once information about n is obtained, a correlation 
with T is possible. In this case, the combination of both parameters can hint towards our 
MSSM RGE scenario. 
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Figure 7. T (top) and R (bottom) as a function of n assuming SM (orange) and MSSM (blue) 
running of mixing parameters. We show the best fit values (left) and full 3(t variation (right), as 
described in the text. 

4.3 Average R and T 

Deep inelastic neutrino-nucleon scattering (DIS), either by charged (CC) or neutral current 
(NC), is the most likely process to take place when an UHE astrophysical neutrino from 
AGN, with energy in the range 10^ < Ey/Ge\ < 10^^, interacts with the nucleons in 
the Earth. As a way to give a more realistic approach to our results we have calculated 
averaged forms of R and T using the DIS detection cross section. 

In order to introduce the averaged forms of the flavour ratios, we have defined an 
averaged transition probability, weighed by the differential CC DIS cross section. This is 
given by the following expression: 

{Pa^{Ql)f = f'dx f\y'^ffl{E,,x,y)P^p{Q^) , (4.18) 

with the integrated CC cross section defined as 



We have used the standard formulae for the differential DIS cross sections [41, 42], with 
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Figure 8. Cross-section-averaged flavour-transition probabilities under MSSM running, eq. (4.18), 
as functions of the threshold transferred momentum, Q^^^. The gray areas are accessible with MSSM 
running when the mixing parameters at production are varied within their current 3(t bounds (see 
text), while in the hatched areas only ^13 and the phases are varied. The bands limited by dashed 
(red) lines are the regions of the ratios accessible with standard oscillations only, i.e., eq. (3.1) with 
Ul = Ui/, when the mixing parameters are allowed to vary in the same manner. In the integral in 
eq. (4.18), we have restricted the transferred momentum values to > Q^^. 



the CTEQ6M parton distribution functions^ (pdf's), fitted within the MS scheme [43]. 
The transferred momentum in a DIS event is = 2mj\fxyE^, with x and y the Bjorken 
scaling parameters^. Its value hes inside the range < < 2mj^E,y, the lower limit 
corresponding to x = y = 0, and the upper limit, to x = y = 1. 

We find that, once we include the pdf's, the average probability favours values of low 
transferred momentum, erasing any effects due to the running. To avoid this, we introduce 
a cut-off in the form of a Heaviside function (Q^ — Q^h) within the integrals in eq. (4.18). 
This artificially restricts the values of to lie above Q^^^, in an attempt to isolate the 
region of parameter space with high values of Q^, where the MSSM-running effects are 
non-negligible. 

In figure 8 we have plotted the cross-section-averaged probabilities, eq. (4.18), as func- 
tions of Q^jj, fixing the neutrino energy at Ei, = 10^^ GeV. As in figure 2, to generate the 

The pdf's are defined for 10"'' < a: < 1 and (1.3)^ < Q^/GeV^ < 10^ outside these ranges, they have 
been extrapolated [44]. 

^The Bjorken scaling parameters are defined as a; = ^ and y = ^^.^ , with pjv and mjv the nucleon 
momentum and mass, respectively. 
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gray regions, we have allowed the mixing parameters at production to vary within their 3a 
experimental bounds, while the CP and Majorana phases have been varied between and 
27r. The hatched regions were generated by varying only ^13 and the phases, while keeping 
all of the other mixing parameters at production at their best-fit values. The dashed (red) 
lines are included for comparison to the standard-oscillation probabilities. In the absence 
of MSSM running effects, the probabilities would have to be confined between these lines. 

From the definition in eq. (4.18), we see that the average probability receives, at a 
given value of Q^^, the accumulated effects from Q'^^ up to (2 x 10^^) m^v GeV^. Thus, if 
we assumed naively that the values of the weights in the integrand, i.e., of the differential 
cross section, remain the same throughout the whole range of integration, then deviations of 
{Pa/S {Qth})'^ from the no-running expectations would occur for all values of Q^j^. However, 
in figure 8, the starting point of the deviations is around 10^ GeV^. For values of Q'^^i ^ 
GeV^, the accumulated effects are not visible, and hence the probabilities under MSSM 
running coincide with the corresponding ones under no running. The reason is that the 
pdf's decrease with growing x, which means that the weights in eq. (4.18) become irrelevant 
at high values of compared to those at low Q^. In contrast, for Q'^^ > 10^ GeV^, the 
weights evaluated at Q'^^ are comparable to the ones evaluated at slightly higher 

values of Q^, with the result that the average probabilities start deviating. Finally, note 
that, with the exception of the differences in the starting points of the deviations, the 
probability regions in figures 2 and 8 have similar shapes and reach the same minimum 
and maximum values. 

The cross-section-averaged flavour fluxes at Earth for the four scenarios, {Ey)), 

are given by eqs. (4.8)-(4.11) after changing — )• (Pa/j {Eu)Y . These flavour fluxes can 
be used to deflne the cross-section-averaged flavour ratios (T) and respectively, by 
eqs. (4.12) and (4.13), with the replacement ^>,_, — >■ ). 

As expected, we find no difference between the cases of no running and MSSM running 
for either (i?) or (T) if a low cut-off or no cut-off at all is imposed. Therefore, in order 
to bring out any differences between the two cases, we show in figure 9 the scatter plot 
of {R) vs. (T), for a fixed energy of = 10^ GeV and a cut-off Ql^ = W GeV^. Since 
this energy is close to the maximum neutrino energy expected, any observed deviations 
of the probabilities will be close to the largest possibly achievable. In accordance to the 
results found in Section 4, the largest potential deviations occur for scenarios A and B, 
corresponding, respectively, to pure electron-flavoured and pure muon-flavoured fluxes. In 
these cases, the stronger deviations are observable in {R). 

4.4 Experimental perspective 

At first sight, the results presented in figure 9 encourage us to look for the MSSM running 
effects in high-energy astrophysical neutrinos within an experimental context. However, 
the requirement of using only scattering events with > Q'^^i — 10'' GeV^ reduces dramat- 
ically the size of the sample of scattering events that could potentially reveal the presence 
of MSSM running in a detection experiment. The situation turns out to be even worse 
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<T> <T> 

Figure 9. Scatter plots of (R) = (l>,,_,+pJ/($^^+p^) vs. (T) = for the four different 

production scenarios A-D (see table 1). Orange points were calculated using standard oscillations 
only, while blue points were calculated under the MSSM running, with an energy = 10^ GeV 
and a cut-off of = 10^ GeV^. 



Scenario 


WB (xlO^) 


BB (xlO^) 


Std. osc. 


MSSM 


Std. osc. 


MSSM 


A 


0.95 - 2.44 


1.12-2.70 


1.16-2.96 


1.34-3.25 


B 


2.51-3.80 


2.22-3.74 


3.70-5.60 


3.32 - 5.51 


C 


1.75-3.16 


1.69-3.27 


2.58-4.66 


2.50 - 4.80 


D 


2.01-3.38 


1.87-3.43 


2.96-4.98 


2.77-5.03 



Table 2. Ranges of values of the integrated number of (i/^ + t^^) expected in IceCube-80 (effective 
volume Voff = 1 km^) in the range 10^ < E^/GcY < 10^^ after t = 15 years of exposure, using 
standard (no-running) and MSSM flavour transitions, for the WB and BB neutrino flux models 
described in Appendix B. A cut-off of Q^j^ = 10^ GeV^ has been applied. The ranges were found 
by varying the mixing parameters at production within their 3a bounds and the phases between 
and 27r. 

if we take into account that at the energy that we have considered, Ei, = 10^ GeV, the 
astrophysical neutrino flux is expected to be very low. 

In order to illustrate this situation, we have calculated the minimum and maximum 
integrated number of muon-neutrinos and anti-neutrinos between 10^ and 10^^ GeV, after 
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15 years of exposure at IceCube, under two assumptions on the UHE diffuse neutrino flux 
from AGN, one by Waxman-Bahcall (WB) and tlie other by Becker-Biermann (BB) (see 
Appendix B), and applying a cut-off of Q"^^ = 10^ GeV^. The minimum and maximum 
numbers of events were obtained by varying the mixing parameters at production within 
their 3a bounds and ah the phases between and 27r. These numbers are displayed in 
table 2. There is a strong trade-off between the value of Q^^ and the expected number of 
neutrino events that exhibit a reconstructed equal to or higher than our threshold: the 
cut-off we have imposed (Q^h ~ GeV^) dramatically decreases the number of events to 
10~^ — 10~^ in 15 years, rendering the potential observation of SUSY running through the 
flavour fluxes null. Thus, detection of a few events would require a Cerenkov detector with 
an unrealistic effective volume roughly 10^ (for the WB flux) to lO'^ (for the BB flux) times 
larger than IceCube's. Using the lowest possible SUSY-enhancing cut-off, = 10^ GeV^, 
would increase the event yield in about three orders of magnitude, but still not enough 
for actual detection. Under these prospects, it makes no sense to calculate the expected 
reconstructed values of the ratios R and T. Therefore, it is evident that IceCube has no 
actual sensitivity to MSSM running in the flavour ratios: we are hampered both by the 
cut-off and by the low flux of UHE neutrinos. 

Although alternative neutrino detection techniques could be considered, such as ra- 
dio [45] or sound detection [46], these would still be limited by the low neutrino flux. 
Furthermore, they would also have to be able to reconstruct the transferred momentum in 
each interaction and the high cut-off needed would also reduce the sample size dras- 
tically. It seems that presently it is technically unfeasible to detect the effects of MSSM 
running in the UHE astrophysical neutrino flavour ratios. 

5 Summary and conclusions 

We have analysed the possibility of observing the renormalisation group running of the 
neutrino mixing parameters in the MSSM. To this end, we have studied the possibility 
of observing distortions in the extragalactic high-energy neutrino flavour ratios, where 
the running effects participate through neutrino flavour-transition probabilities. With the 
inclusion of SUSY effects, modifications to these probabilities with respect to their standard 
values grow with the transferred momentum in neutrino-nucleon deep inelastic scattering, 
Q^, which is largest for UHE neutrinos from AGN. 

It has been observed that the transition probability under MSSM running starts to 
differ from the one with no running at = 10^ GeV^, which corresponds to our choice 
for the SUSY scale Asusy- The maximum variation is in the order of 50% for P^^-- One 
important feature of these SUSY effects is that they are mainly controlled by the Dirac 
and Majorana CP phases, and by ^13. We have found that in order to enhance the SUSY 
effects two conditions must be satisfied: a large value of tan/3 and a large value for the 
sum of the neutrino masses. We have set tan/3 = 50 and Tiin,^. < 1.3 eV, following the 
latest WMAP-only bound. 

We have studied R = $^,+p,/^>^,+p, and T = ^u^+pj (^>^,+p, + ^>^^+p^ + ^u^+p^), 
taking into account the MSSM effects on the oscillation probabilities, for four different pro- 
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duction model scenarios, ($0^+^^ : $°^+p^, : ^°,+p,) = (1 : : 0), (0 : 1 : 0), (1/2 : 1/2 : 0), 
and (1/3 : 2/3 : 0), considering the full variation of the neutrino mixing parameters within 
their current 3a allowed ranges. We have found that, under MSSM running, for ~ 10^^ 
GeV^, or even for smaller values, R is able to reach, within the (1:0:0) production sce- 
nario, values 18 times larger than the corresponding ones in the no-running case, while for 
T the differences are not nearly as significant. In the remaining neutrino production scenar- 
ios, the differences between the MSSM running and no-running cases are less pronounced, 
more so for (1/2 : 1/2 : 0) and (1/3 : 2/3 : 0). 

In order to give a more realistic view of our observations on R and T, we have built 
average forms of these observables using as averaging weights the DIS cross section (i.e., 
the cross section at detection). We have noted that, in order to obtain differences between 
the MSSM-running from the standard behaviour in our predictions of (R) and (T), it is 
necessary to impose a cut-off condition of > Q^j^ = 10''' GeV^. Otherwise, if we consider 
lower values of the threshold Q"^^, the differences completely disappear. This result is 
discouraging for going forward in the search for MSSM effects in the flavour ratios, since 
two factors suppress the event rate at a large neutrino telescope: the expected low UHE 
astrophysical neutrino flux and the fact that the event sample at high Q"^^ is too low. We 
have estimated the number of events at IceCube to be 10~^ — 10~^ after 15 years, for a 
cut-off of Qt\ = 10^ GeV2. 

Therefore, although there are corrections to the flavour transition probability due to 
the MSSM running of the mixing parameters, we are compelled to conclude that it is not 
feasible, with the available and foreseeable technology, to detect the presence of such effects 
in the flavour ratios of UHE astrophysical neutrinos. 
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A Neutrino flavour-transition probability with diff"erent mixing matrices 

The probability of a neutrino produced with flavour a interacting with flavour /3 can be 
expressed as: 

3 

Pa(5 = ^9aig'i3i , (A.l) 
1=1 

where the coefficients gai depend only on (6*12,^13,6*23,(5) and the g'j^- depend only on 
(Ci2) Ci3 5 C23 5 (^i) and, through the RGEs, on Q. Their expressions are shown in table 3. 
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Coefficient name 


Expression 


Best-fit value 


Qel 


„2 „2 


U.O ( 


9e2 


2 2 


(J. 32 


5e3 


2 


0.01 




4^ / O O O X^ 
/ (.(712, c/i3, ^23, Oj 


O.i r + 0.045C5 


9n2 


/(012+7r/2,%,^23,5) 


0.37 - 0.045C5 




„2 2 


0.46 


5t1 


/(ei2+vr,ei3,^23+vr/2,5) 


0.15 - 0.045C5 


5t2 


/ (012 + 7r/2, 013, 023+^/2, 5) 


0.32 + 0.045C5 


5r3 


6*13 023 


0.53 



Table 3. g-functions used to calculate the flavour-transition probability, eq. (A.l). The ex- 
pressions for arc obtained from the corresponding expressions for g^i shown here, after the 
replacements 9ij — >■ Qj and (5 — >■ The auxiliary function / is defined as f {6i, 02,03,6) = 



Co„ So 



Thus, we see that the probability can be written as a sum of terms, each of which is 
the product of a function of only standard mixing parameters times a function of only 
renormalized mixing parameters. 

Note that the standard probabilities can be recovered simply by taking — )• Oij and 
(5i 5, which is equivalent to the replacement g'^^ gpi in eq. (A.l). 

The standard coefficients [gai — gaj) shown in eqs. (3.2)-(3.4) can then be written as: 



9el 


- 9e2 


9el 


- 9e3 


ge3 


- 9e2 




- g^l2 


g^ll 


- gtJ.3 


9tJ.3 


- g^l2 


9t1 


- gr2 


9t1 


- grs 


9t3 


- gr2 



2 2 



C0i3C2ei2 

C2ei3 
-C2013 +^^12^013 

C5S2ei2 52023 S9l3 - C2012 (C2e23 + ^lia^ln) 

1 

2' 

-C<5S26»i2'S26»23'5ei3 + C26»i2 (£26*23 



T^C5S2ei2 52^23 ■5613 - cL^L + (c2ei3 + cL^D ■'^82 

2 2 \ 
%3%3/' 



1 



C5S26»i2S26»23 56*13 + ^923^91 



(c2ei3 +cLsL)ci 



1 



:C5S29i2 52623 •Sfii 



2 2 / 2 2 \ 2 

'^023'^6'l2 + l'=26»i3 + ^012^ en J '^023 



(A.2) 
(A.3) 
(A.4) 
(A.5) 

(A.6) 

(A.7) 
(A.8) 
(A.9) 

(A.IO) 



B Number of neutrinos expected at IceCube 

Adapting the expressions from [47], the number of CC interactions initiated by a down- 
going'^ astrophysical flux of a-flavoured (anti-)neutrinos in a detector of effective volume 

^ Since at energies E„ > 10^ GeV energy losses of Earth-traversing neutrinos in NC VN interactions 
become important, i.e., the shadow factor [42] becomes < 0.1 and the Earth turns opaque for upgoing 
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Ves = 1 km'^ and opening solid angle Q = 5.736 sr (corresponding to a zenith angle of up 
to 85°) can be calculated as 

AtCC = turVesn / ' dE, 4c (E,) (|.,_, (^.))cc ^-..n (^v) , (B.l) 

where E"™'" = 10^ GeV and E^^^^ = 10^^ GeV are, respectively, the minimum and max- 
imum neutrino energies that we have considered, t = 15 years is the exposure time, 
riT = 5.1557 X 10^^ cm~^ is the number density of targets (nucleons) in ice, and a^Q 
is the (anti-) neutrino CC scattering cross section off an isoscalar target, eq. (4.19). 

The flux ^u^ii is the all-flavour diffuse flux of astrophysical neutrinos, for which we 
have used two models: a conservative one, by Waxman and Bahcall (WB) [48] and an 
optimistic one, by Becker and Biermann^ (BB) [49]: 

$WB ^ (^^/GeV)-2 GeV"^ cm'^ s'^ si'^ (B.2) 
(Ei.) =^ 5 X 10^2 (E^/GeV)-'^-^ GeV"^ cm'^ s'^ sr^^ . (B.3) 

Both fluxes lie below the upper bound curves reported by different neutrino detection 
experiments: AMANDA-II [52, 53], ANITA-II [54], Auger [55], RICE [56], IceCube 22- 
strings [52], and a preliminary IceCube 40-strings bound [57, 58]. In particular, the BB 
model parameters were chosen to yield the maximum allowed flux within this model^. 
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